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Abstract 

We investigate theoretically baryon systems made of three hadrons which contain one nucleon and one 
D meson, and in addition another meson, D,K or K. The systems are studied using the Fixed Center 
Approximation to the Faddeev equations. The study is made assuming scattering of a if or a J? on a DN 
cluster, which is known to generate the A c (2595), or the scattering of a nucleon on the DD cluster, which 
has been shown to generate a hidden charm resonance named X(3700). We also investigate the configuration 
of scattering of N on the KD cluster, which is known to generate the D* (2317). In all cases we find bound 
states, with the NDK system, of exotic nature, more bound than the KDN . 

PACS numbers: 14.40.Lb; 14.20.Lq; 21.45.-v 
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I. INTRODUCTION 



While the three baryon system has been a subject of intense theoretical study [H-El], ^ has 
only been recently that attention was brought to systems with two mesons and one baryon, with 
unexpected results. Indeed, states with two pseudoscalar mesons and one baryon were studied in 
0], combining Faddeev equations and chiral dynamics, by means of which the low lying excited 
J p = 1/2+ A and S states of the Particle Data Book (PDG) [5|] were described. The same happened 
with the low lying J p = l/2 + N* states 0]. Independently, and using variational techniques, a 
N* state around 1920 MeV was predicted in 0] as a molecule of NKK, corroborated in coupled 
channels Faddeev equations in ja] and js)]. A different calculation also predicts a quasibound 
irKN system, with the difference that the Ntt interaction is in p-wave [lOj]. Systems of bound 
or quasibound three mesons did not wait long and in lj| the X(2175) (now the 0(2170)) was 
explained as a resonant <f>KK system, also described as such in a phenomenological way in [12|. In 
a similar way, the -fT(1460) is explained as a KKK state in (l3| and more recently the 7r(1300) is 



described as a irKK molecule in 14] 



The charm sector has not yet been explored for such three body systems and this is the first 
incursion in that field. For this purpose we have selected systems that have a nucleon and a 
D meson. The DN system, in collaboration with coupled channels, leads to the formation of a 
dynamically generated state, the A c (2595) [l5l4l8|. On top of it we add a K, K or D meson and we 
study the stability of the s yste m. The case of scattering of N on the DK cluster, which is known 
to generate the £>* (2317) )±dt[~ 



a bound state in isospin 1=0 



4211] is also considered. On the other hand, the DD system leads to 
211 ]. which might have already been observed 22|, in analogy to the 
/o(980) which is largely made of KK [2314291] . We add a nucleon to it and study the interaction 
of the three body system. The system obtained would be the analogous to the one found in 0-0] 
as a NKK. In all cases we find bound or quasibound states with masses around 3100 MeV in the 
first cases and 4400 MeV in the case of NDD. 

The calculations are made using the Fixed Center Approximation (FCA) to the Faddeev equa- 
tions. The method has proved to be rather reliable for cases like -fT-deuteron scattering very close 
to threshold (see [13] and the earlier work [H|]). I n a closer problem to the present one, diverting a 
bit from threshold in the bound region, the FCA has been applied to the study of the NKK system 
in [32| and the results compare favorably with those of the Faddeev approach in [s| and those of 
the variational approach in M|. Similarly the FCA has been applied to an analogous problem, the 
one of the KNN system [33j and a detailed exposition has been made in the Introduction of the 
approximations involved in ordinary Faddeev calculations, which can induce uncertainties as large 
as those of the FCA for different reasons. Yet, it has been reassuming to see in that paper that 
the results of the FCA are qualitatively in agreement with those of other app roaches 3414391] , and 
remarkably similar to those obtained in the variational approach of 37|, [38J] which use the same 



input as in [33j. The differences between [37|, |38J] and [33J] are at the level of a few MeV in the 



binding, an accuracy which is far more that sufficient in the exploratory work that we carry on the 
systems under consideration, which are studied here for the first time. 



II. MULTI-BODY INTERACTION FORMALISM 



The Faddeev equations under the FCA are an effective tool to deal with multi-hadron interaction 
They are particularly suited to study system in which a pair of particles cluster 
together and the cluster is not much modified by the third particle. Even if there is a sizeable 
modification of the cluster, the method is useful in combination with information from other sources 



on how the cluster can be modified by the presence of the third particle 33)]. The assumption of a 
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small modification of the cluster wave function seems to imply that the mass of the third particle 
should be smaller than that of the cluster components, but this can also happen if the cluster is 
strongly bound, independent of the masses. In any case, in an exploratory study of the systems 
under consideration, where uncertainties of even 50 MeV are readily acceptable, the FCA proves 
to be a sufficiently accurate tool, and relatively simple to use, once comparison with more accurate 
tools has shown that the same results are obtained within a few MeV of difference. 

Following 32|, 33|, 40-42], we will apply the FCA to study the charm sector. The FCA approxi- 



mation to Faddeev equations assumes a pair of particles (1 and 2) forming a cluster. Then particle 
3 interacts with the components of the cluster, undergoing all possible multiple scattering with 
those components. This is depicted in Fig. [TJ With this meaning of the FCA, it is easy to write 
the equations. For this one defines two partition functions Ti, T 2 which sum all diagrams of the 
series of Fig. Q] which begin with the interaction of particle 3 with particle 1 of the cluster (Ti), or 
with the particle 2 (T2). The equations then read 

Ti = t 1 +t 1 G T 2 , (1) 
T 2 = t 2 + t 2 G Q T l} (2) 
T = Tx+T 2 , (3) 

where T is the total three-body scattering amplitude that we are looking for. The amplitudes t\ 
and t 2 represent the unitary scattering amplitudes with coupled channels for the interactions of 
particle 3 with particle 1 and 2, respectively. Besides, Go is the propagator of particle 3 between 
the components of the two-body system, which will be discussed below. 

Now we turn to the amplitudes corresponding to single-scattering contribution. One must take 
into account the isospin structure of the cluster and write the t\ and t 2 amplitudes in terms of the 
isospin amplitudes of the (3,1) and (3,2) systems. Details can be seen in [32j, |4l|. In the present 
case this is particularly easy. Whether we have the KDN, KDN, NDK or NDD system, where 
the first particle is labelled 3 and the last two particles are making the cluster (particles 1, 2), the 
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bound state of 1, 2 is in I = from former studies and the total spin is then I = ^. Then, for all 
four cases we find 
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in which t^ 1 is the two-body unitary scattering amplitude of isospin 1 = 1 between particle 3 
and 1 evaluated with its coupled channels, and similarly for the other cases. We show below the 
explicit evaluation for the NDD case. The derivation for the other systems is identical. Here we 

I ndd = 1/2- We have 



take the case of I D d = and Itotai 
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where the notation of the states followed in the terms is |(1, 1), — I >= \(I 31 , 1 31 ), 1^ > for £31, and 
|(/ 32 ,/f),/i>fort 3 2. 

Because we follow the normalization of Mandl and Shaw [43] which has different weight factors 
for the meson and baryon fields, we must take into account how these factors appear in the single 
scattering and double scattering and in the total amplitude. This is easy and is done in detail 
in (32I . 41 ] for the two different cases, when the cluster is a baryon (KDN, KDN) or a meson 
(NDD, NDK). We show below the details for the case of a baryonic cluster and a meson as 
scattering particle. Let us assume that particle 1 of the cluster is a meson and particle 2 is a 
baryon. Particle 3 is also a meson. In this case, following the field normalization of 43J we find 



4 



for the S matrix of single scattering, 
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1 1 



1 



2M 2 2M 2 



v 2 V2^y^J 3 y 2E 2 i 

x (27r) 4 <5(fc + k R - k' - k' R ), 



2E' 2 



(8) 



(9) 



where, fc, A;' fcp) refer to the momentum of initial, final scattering particle (R for the cluster), 
V is the volume of the box where the states are normalized to unity and the subscripts 1, 2 refer 
to scattering with particle 1 or 2 of the cluster. 

The double scattering diagram, Fig. CD (b), is given by 



S {2) = - i(27r) 4 5(k + k R -k' - k' R ) 



1 1 
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1 



(2tt) : 
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2M 2 \2M 2 
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(10) 



mi + i e 



where F R (q) is the cluster form factor that we shall discuss below. 

Similarly the full S matrix for scattering of particle 3 with the cluster will be given by 

S = - i T (2n) A 5(k + k R — k' — k' R )^ 



1 



1 2M R 2M R 



(11) 



2o; 3 J2U, V 2E R V 2E[ 



R 



In view of the different normalization of these terms, we can introduce suitable factors in the 
elementary amplitudes, 



1 



2mi 



ti, t 2 — t 2 , 



where we have taken the approximations, suitable for bound states, 
sum all the diagrams by means of 

ti + t 2 + 2 ti t 2 Gq 
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(12) 
1, and 

(13) 



The function Go in Eq. f|13f) is given by 



G (s) 



d 3 q 
(2tt) 3 
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(14) 



where F R (q) is the form factor of the cluster of particles 1 and 2. We must use the form factor of 
the cluster consistently with the theory used to generate the cluster as a dynamically generated 
resonance. This requires to extend into wave functions the formalism of the chiral unitary approach 
developed for scattering amplitudes 44-53j]. This work has been done recently in 54, 55] and the 
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expression for the form 
reproduce below 



factors are given in section 4 of 



551 ]. which we use in the present work and 
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(16) 



where E\ and E 2 are the energies of the particles 1, 2 and Mr the mass of the cluster. The 
parameter A is a cut off that regularizes the integral of Eqs. (|15j) and (I16j) . This cut off is the same 
one needed in the regularization of the loop function of the two particle propagators in the study 
of the interaction of the two particles of the cluster 55J] . We take these values of A such as to get 
the bound states of A c (2595) from 0, L>* (2317) from and AT(3700) from 0. 

When the cluster is a meson and the scattering particle a baryon, the solution is given by the 
same Eq. (I13|) . but now 
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In addition, q°, the energy carried by particle 3 in the rest frame of the three particle system, is 
given by 



9° 00 



s + mi 



Ml 



2^ 



(19) 



Note also that the arguments of the amplitudes Tj(s) and tj(sj) are different, where s is the 
total invariant mass of the three-body system, and Sj are the invariant masses in the two-body 
systems. The value of Si is given by [41] 



m i + m i + 



(Ml + mj 



m|)(s 
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Ml) 



2Ml 



= 1,2, i^j) 



(20) 



where mi, (I = 1,2,3) are the masses of the corresponding particles in the three-body system and 
Mr the mass of two body resonance or bound state (cluster) . 



III. THE CASE OF KDN INTERACTION 

Our strategy proceeds as follows: first we generate the resonance or bound state in the com- 
pound system and determine the value of the parameter A, then calculate the form factor and Go 
propagator and take the t\ and t 2 amplitudes from the unitary coupled channel approach, finally 
the total scattering amplitude T is evaluated. For the KDN scattering, the first thing we do is 
to reproduce the work of in coupled channels for the DN system. The coupled channels 

used are 7rE c , DN, r]A c , K3 c , Kr! c , D S A, 7]'A C and the dynamics is obtained from the exchange 
of vector mesons between the pseudoscalar meson and the baryon. This procedure, based on the 



local hidden gauge approach [57H59J, leads to the chiral Lagrangians in the SU(3) sector. One 



gets as dynamically generated resonance the A c (2595), which couples most strongly to the DN 



6 




2600 2700 

s 1/2 [MeV] 



3000 



FIG. 2: Imaginary part of the DN amplitude for isospin 1 = 0. 



channel and is interpreted as a DN bound state. As shown in Fig. El the method generates a 
pole in the DN scatte ring amplitude in I = in the first Riemann Sheet at (2595 + iO) MeV. 
Since the works of use dimensional regularization for the loops, and we need a cut off to 

obtain the wave function and form factor, the equivalent cut off must be obtained. There are three 
methods to do this. One of them is to compare the value of the G propagator (loop function of 
two particles propagator which appears in the Bethe-Salpter equation T = V + VGT) at threshold 
using the dimensional regularization formula 48] with the one of the cut off which can be taken 



from [45J ( [24 ] for meson- meson interaction) or the analytic expression in [60]. Another method is 



to compare the pole position using different regularization schemes. The third one is to use the 
relation between the subtraction constant a(fj,) and the cut off A of Eq. (52) in 61]. The best 
fitting results by these methods give us a value of A = 973 MeV. The A c (2595) form factor using 
this cut off is shown in Fig. [3j In the next step we evaluate Go by means of Eq. (I14D and we show 
its real and imaginary parts in Fig. 01 

In the final calculation we also need to know the two-body unitary scattering amplitudes in 
different isospin states. For the KDN interaction, the amplitudes of KD and KN are needed. For 
the KD interaction we have taken the results of [2l| . On the other hand, the KN scattering has 
been evaluated using the chiral unitary approach of [45] with the dimensional regularization scheme 
of [!§], using n = 630 MeV, aj(/x) = —1.84 for all channels. This scheme leads to the generation of 
the two A(1405) states reported in [5ol ]. 

In Fig. [5]we show the results of |T| 2 for the -KA C (2595) scattering. We find a peak around 
3150 MeV, slightly above the threshold of the A c (2595) + K mass (3088 MeV) and below the 
threshold of the KDN system (3298 MeV). The width of the peak is about 50 MeV, which 
indicates the width of the state that we obtain. In our study of the system KDN, where we have 
chosen the DN system in / = since this is the channel with strong attraction leading to the 
A c (2595) resonance, the quantum numbers of the KDN state are C = +1,5" = —1 and J p = | + 
since we only consider the interaction among the components in L = 0. The mass of this state 
is very close to that of the 3 C (3123), of unknown spin and parity, but its decay width is larger 
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FIG. 3: Form factors of A c (2595) and L»* (2317). 




s 1 ^ [MeV] 

FIG. 4: Real(solid line) and imaginary(dashed line) parts of Go in KDN. 

than that of the S c (3123), 4 ± 4 MeV The larger width, tied to the vrS decay of the KN 
system, seems unavoidable, and this indicates that the resonance that we find could most probably 
correspond to a resonance not yet found. 
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IV. INVESTIGATING THE NDK INTERACTION 



The two body KD and DN interactions were studied by the coupled channel Bethe-Salpeter 
equations in 

[El and 

respectively. It was found that the resonance 7}* (2317) is dynamically 



generated in 7 = from KD scattering and the A c (2595) is produced in 7 = from the DN 
interaction, as we mentioned above. Hence there are two possible cases of three body scattering in 
the NDK system. One is the N — (DK)jj* o (^sn) and the other one is the K — (7JW)a c ( 2 595)- 

First, we are going to investigate the three body scattering for the N — (DK)n* (9.^7) ■ In 
order to calculate this, in a first step we obtain the DN amplitude, t\ from (la . l56| and the KN 



amplitude, ty, from the chiral unitary approach of |45l |. For the DN matrix element the result of 
[l6| is reproduced and the imaginary part of the DN matrix element that we obtain is shown in 
Fig. [21 In the case of the KN system, the interaction is repulsive in 7 = 1 and has vanishing 
interaction in 7 = 0. We take the parameters for the loop function from (48| . For the next step, 
in order to get the total scattering amplitude T, we need to know the form factor for 7)* (2317) 
and Gq(s). Here the cut off is determined by comparing the value of the G function that one 



obtains from the dimensional regularization [48[ and the cut off scheme 45j at threshold. In this 



way A = 900 MeV is obtained for the cut off of the 7)* (2317) form factor. The form factor of the 
7)*q(2317) is plotted in Fig. [3] and the Gq(s) function for this case is represented in Fig. El 

Using the aforementioned total three-body scattering amplitude T, we obtain |T| 2 for the 
V7?* (2317) scattering shown in Fig. [71 We found a peak around 3050 MeV which is about 
200 MeV below the 7?* (2317) and N threshold. This reflects the strong attraction in the DN 
system that leads to the A c (2595). The width of the state is smaller than 10 MeV. We do not find 
a counterpart in the PDG and the quantum numbers, with positive strangeness, correspond to an 
exotic state. 

As an alternative possibility, the three body scattering in the KDN system can also be formed 
as K — (TW) Ac (2595)j with the A c (2595) assumed as a two body cluster. Now we need the two body 
matrix elements KD and KN. The latter one is calculated using the chiral unitary approach that 
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FIG. 6: Real(solid line) and imaginary(dashed line) parts of Go in N — DK . 
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FIG. 7: Modulus squared of the iVD* (2317) scattering amplitude. 



was mentioned before. The KD matrix element, as also mentioned above, is investigated in [2l|. In 
order to calculate the form factor for the A c (2595), the cut off is needed. With the aforementioned 
strategy the value of A = 973 MeV is obtained and it was already used in section IIII1 Using this 
cut off value, the Gq{s) function is the same as in section [nil when we had scattering of a K on 
the cluster of A c (2595), since the masses of K and K are the same. The Gq(s) function is, thus 
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FIG. 8: Modulus squared of the ifA c (2595) scattering amplitude. 

the one plotted in Fig. [H Ultimately the total three-body scattering amplitude for the if A c (2595) 
scattering is evaluated and the result of \T\ 2 is represented in Fig. [8j There is an explicit narrow 
peak at 3100 MeV. 

One may try to investigate the structure of the peak in Fig. [8j but it would only distract us 
from the main point, which is that the weight of \T\ 2 in Fig. [8jis very small compared with the 
one in Fig. [7] for the N — (DK)d* (2317) configuration. Note that a proper comparison requires to 
take into account the different field normalizations. Indeed, the S matrix for N — (-D-^Qd* (2317) 
goes as 



S ~ 1 - %T 



V 2a ^ (2317) \l 2j D lo (2Zn) 



while for K — (£W)a c (2595) it goes as 



2M N 2M N 
2En \ 2E* 



N 



(21) 



S ~ 1 - %T 



1 2M Ac 2M Ac 



^2^ yf2J^ y 2E Ac V 2E' A 



(22) 



Hence, the proper comparison is 



2m T 



in the first case versus -^-^ in the second, or T(K(DN)) 



versus 



x T{N{DK)). Considering this, the strength of the peak for \T{K(DN))\ 2 is about a 

factor 90 smaller than for \T(N(DK))\ 2 . This means that the K(DN) configuration in the wave 
function of the KDN system has a very small weight. Hence, we predict a bound state of NDK 
mostly made of a N orbiting around a bound DK cluster forming the -D* (2317). 
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FIG. 9: Real(solid line) and imaginary(dashed line) parts of Go in NDD. 



V. NDD INTERACTION RESULTS 



USB and a resonance called X(3700), 



The two-body DD interaction was investigated in 
was dynamically generated. This resonance would be the analogue one to the /o(980) which is 
basically a KK bound state (23l-[29|. In our procedure we also reproduce this DD state successfully, 
getting the pole as (3718 + iO) MeV with the same parameters as in 21]. We take a value of 
A = 850 MeV from 0, which is consistent with the methods mentioned above. Then we can 
calculate the form factor of the X(3700) with Eq. (|15p by means of this cut off. Next we evaluate 
the Go by means of Eq. (|18p . for N propagating between the D and D, and its results are shown 
in Fig. El 

The nucleon interaction with the D, D mesons was studied by the coupled channels two-body 
scattering equations in 16L 17]. For the DN scattering amplitude, as mentioned before, we followed 



the procedure of Iff]. For the 1 = 1 sector there are eight coupled channels and we have used the 
same parameters as in / = which reproduced the A c (2595) resonance. The DN interaction, 



which is similar to the KN channel [16j j. is repulsive in I = 1 and vanishes for 1 = 0. As in 171 ]. 
we took the same parameter as for the DN interaction. Finally we obtain the T matrix, for the 
NDD interaction by means Eq. (|3|), and show the results of |T| 2 in Fig. [TPl From this figure we 
can see that there is a clear peak of |T| 2 around 4400 MeV and the width is very small, less than 
10 MeV. The peak appears below the NDD and iV X (3700) thresholds and corresponds to a bound 
state of iVA(3700). This would be a hidden charm baryon state of J p = \ + which appears in the 
same region of ener gies another hidden charm states of J p = h obtained from the DA r , Z)£ r 
and D*\ 



D* Yj c in [63] , 1 641 . In these latter works some reactions were suggested to observe those 



states in future Facilities. The same or similar reactions could be used to observe these states of 
positive parity. 
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FIG. 10: Modulus squared of the NX (3700) scattering amplitude. 
VI. CONCLUSIONS 

We have investigated three body systems that have one D meson or DD, together with one 
baryon. The systems are KDN, NDK and NDD. We have investigated the interaction of the 
systems taking into account that a couple of the particles in each case are strongly bound, forming 
dynamically generated states. We let then the third particle interact with the components of 
this cluster. Concretely, for KDN we study the scattering of K with the D and N components 
of the cluster of DN that makes the A c (2595) resonance. In the second case, NDK, we find 
that the important configuration corresponds to N scattering over the cluster of KD that makes 
the D*q(2317). In the case of the NDD we look at N scattering on the DD cluster that is 
supposed to generate a bound state called X(3700), in analogy to the /o(980) which is mostly 
a bound state of KK. In all cases we find bound or quasibound states, relatively narrow, with 
energies 3150 MeV, 3050 MeV and 4400 MeV, respectively. All these states have J p = 1/2+ and 
isospin 1 = 1/2 and differ by their charm or strangeness content, S = — 1,C = 1, S = 1,C = 1, 
S = 0,C = 0, respectively. The first state could perhaps be associated to the S(3123), which has 
unknown J p , but the width obtained is a bit too large. The second state is of exotic nature and 
there is no counterpart in the PDG. The third state is a regular N* state as to quantum numbers, 
but it contains hidden charm. It lies in an energy region where baryon states have not yet been 
investigated. We are thus making predictions in the frontier of this field and hope that with the 
coming Facilities of FAIR, or the BELLE upgrade, such states can be systematically studied and 
that our predictions can be confirmed. 
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